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ABSTRACT 



We report the first high spatial resolution submillimeter continuum observations of the Sagittarius B2 cloud complex using the 
Submillimeter Array (SMA). With the subarcsecond resolution provided by the SMA, the two massive star-forming clumps Sgr B2(N) 
and Sgr B2(M) are resolved into multiple compact sources. In total, twelve submillimeter cores are identified in the Sgr B2(M) re- 
gion, while only two components are observed in the Sgr B2(N) clump. The gas mass and column density are estimated from the dust 
continuum emission. We find that most of the cores have gas masses in excess of 100 Mq and column densities above 10'^ cm"'. The 
very fragmented appearance of Sgr B2(M), in contrast to the monolithic structure of Sgr B2 (N), suggests that the former is more 
evolved. The density profile of the Sgr B2(N)-SMA1 core is well fitted by a Plummer density distribution. This would lead one to 
believe that in the evolutionary sequence of the Sgr B2 cloud complex, a massive star forms first in an homogeneous core, and the rest 
of the cluster forms subsequently in the then fragmenting structure. 

Key words. ISM: clouds: radio continuum — ISM: individual objects (Sgr B2) — stars: formation 



1. Introduction 

The Sagittarius B2 star-forming region is located ~ 100 pc from 
Sgr A*, within the ~ 400 pc wide dense central molecular zone 
(CMZ) of the Galactic center, at a distance of ~8 kpc from the 
Sun (Reid et al. 2009). It is the strongest submillimeter contin- 
uum source in the CMZ (Schuller et al. 2009). It contains dense 
cores, Sgr B2(N) and Sgr B2(M), hosting clusters of compact Hii 
regions (Gaume et al. 1995; de Pree, Goss & Gaume 1998). It has 
been suggested that these two hot cores are at different evolution- 
ary stages (Reid et al 2009; Lis et al. 1993; Hollis et al. 2003; Qin 
et al. 2008). Spectral observations in centimeter and millimeter 
regimes have been conducted towards Sgr B2 (e.g. Carlstrom & 
Vogel 1989; Mehringer & Menten 1997; Nummelin et al. 1998; 
Liu & Snyder 1999; Hollis et al. 2003; Friedel et al. 2004; Jones 
et al. 2008; Belloche et al. 2008), suggesting that Sgr B2(N) 
is chemically more active. Nearly half of all known interstellar 
molecules were first identified in Sgr B2(N), although sulphur- 
bearing molecules are more abundant in Sgr B2(M) than in 
Sgr B2(N). 

The differences between Sgr B2(N) and Sgr B2(M), in terms 
of both kinematics and chemistry, may originate from differ- 
ent physical conditions and thus different chemical histories, or 
may simply be an evolutionary effect. A clearer understanding 
of the small-scale source structure and the exact origin of the 
molecular line emission is needed to distinguish between these 
two possibilities. In this Letter, we present high spatial resolu- 
tion submillimeter continuum observations of Sgr B2(N) and 
Sgr B2(M), using the SMAQ. The observations presented here 



resolve both of the Sgr B2 clumps into multiple submillimeter 
components. Combining with the SMA spectral line data cubes 
and ongoing Herschel/HIFI complete spectral surveys towards 
Sgr B2(N) and Sgr B2(M) in the HEXOS key project, these ob- 
servations will help us to answer fundamental questions about 
the chemical composition and physical conditions in Sgr B2(N) 
and Sgr B2(M). 

2. Observations 

The SMA observations of Sgr B2 presented here were carried 
out using seven antennas in the compact configuration on 2010 
June 11, and using eight antennas in the very extended con- 
figuration on 2010 July 11. The phase tracking centers were 
^(J 2000.0) = 17^47" 19.883% (J (J 2000.0) = -28°22'18.4" for 
SgrB2(N) and a(J 2000.0) = 17''47'"20.158%(5 (J 2000.0) = 
-28°23'05.0" for Sgr B2(M). Both tracks were observed in 
double-bandwidth mode with a 4 GHz bandwidth in each of the 
lower sideband (LSB) and upper sideband (USB). The spectral 
resolution was 0.8125 MHz per channel, corresponding to a ve- 
locity resolution of ~0.7 km s"'. The observations covered rest 
frequencies from 342.2 to 346.2 GHz (LSB), and from 354.2 to 
358.2 GHz (USB). Observations of QSOs 1733-130 and 1924- 
292 were evenly interleaved with the array pointings toward 
Sgr B2(N) and Sgr B2(M) during the observations in both con- 
figurations, to perform antenna gain calibration. 

For the compact configuration observations, the typical sys- 
tem temperature was 273 K. Mars, QSOs 3c454.3, and 3c279 
were observed for bandpass calibration. The flux calibration was 



' The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
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based on the observations of Neptune (~1.1")- For the very ex- 
tended array observations, the typical system temperature was 
292 K. Both 3c454.3 and 3c279 were used for bandpass cali- 
bration, and Uranus (~1.7") was used for flux calibration. The 
absolute flux scale is estimated to be accurate to within 20%. 

The calibration and imaging were performed in Miriad 
(Sault, Teuben & Wright 1995). We note fliat there are spectral- 
window-based bandpass errors in both ampUtude and phase on 
some baselines in the compact array data, which were corrected 
by use of a bright point source using the BLCAL task. The sys- 
tem temperature measurements for antennas 2 and 7 in the very 
extended array data were not recorded properly and were cor- 
rected using the SMAFIX task. As done by Qin et al. (2008), we 
selected line free channels which may contain some contribu- 
tions from weak and densely spaced blended lines, although at 
least obvious line contributions were excluded. The continuum 
images were constructed from those, combining the LSB and 
USB data of both the compact and very extended array observa- 
tions. We performed a self-calibration on the continuum data us- 
ing the model from 'CLEANed' components for a few iterations 
to remove residual errors. Using the model from 'CLEANed' 
components, the self-calibration in our case did not introduce 
errors in the source structure, but improved the image quality 
and minimized the gain calibration errors. The final images were 
corrected for the primary beam response. The projected base- 
lines ranged from 9 to 80 kA in the compact configuration and 
from 27 to 590 kA in the very extended configuration. The result- 
ing synthesized beam is 0"4x0!'24 (PA=14.4°) using uniform 
weighting, and the Icr rms noise levels are 21 and 31 mJy for 
Sgr B2(M) and Sgr B2(N) images, respectively. The difl'erence 
in rms noise is caused by having more line free channels for con- 
tinuum images in Sgr B2(M) (2417 channels) than in Sgr B2(N) 
(740 channels). Since there are no systematic offsets between 
the submillimeter and cm sources, we beUeve that the absolute 
astrometry is as good as 0.1". 

3. Results 

Continuum images of Sgr B2(N) and Sgr B2(M) at 850 /im 
are shown in Figure 1 . Multiple submillimeter continuum cores 
are clearly detected and resolved towards both Sgr B2(M) and 
Sgr B2(N) with a spatial resolution of 0"4x0!'24. Unhke in ra- 
dio observations at 1.3 cm with a comparable resolution, which 
detected flie UC Hii regions Kl, K2, K3, and K4 (Gaume et al. 
1995), only two submillimeter continuum sources, SMAl and 
SMA2, are observed in Sgr B2(N) (Fig. 1). The bright compo- 
nent Sgr B2(N)-SMA1 is situated close to the UC Hn region 
K2, and Sgr B2(N)-SMA2 is located ~5" north of Sgr B2(N)- 
SMAl. The observations showed that large saturated molecules 
exist only within a small region (< 5") of Sgr B2(N) called 
file Large Molecule Heimat, Sgr B2(N-LMH) (Snyder, Kuan & 
Miao 1994). Our current observations indicate that Sgr B2(N- 
LMH) coincides with Sgr B2(N)-SMA1. Sgr B2(N)-SMA2 was 
also detected in continuum emission at 7 mm and 3 mm (Rolffs 
et al. 201 1; Liu & Snyder 1999) and molecular fines of CH3OH 
and C2H5CN at 7 mm (Mehringer & Menten 1997; HoUis et 
al. 2003). Lower resolution continuum observations at 1.3 mm 
(see note to Table 1 of Qin et al. 2008) suggested that the Sgr 
B2(N)K1-K3 clump could not be fitted with a single Gaussian 
component, and another component existed at Sgr B2(N)-SMA2 
position. Our observations here resolved out Sgr B2(N)-SMA2, 
confirming it to be a high-mass core. 

The continuum image of Sgr B2(M) (Fig. 1, right panel) 
shows a complicated morphology, with a roughly north-south 



extending envelope encompassing several compact components. 
In total, twelve submillimeter sources are resolved in Sgr B2(M). 
Using the Very Large Array (VLA), Gaume et al. (1995) de- 
tected four bright UC Hii regions (F1-F4) at 1.3 cm within 
2" in Sgr B2(M). The highest resolution 7 mm VLA image 
with a position accuracy of 0.1" (de Pree, Goss & Gaume 
1998) resolved nineteen UC Hn regions in the central re- 
gion of Sgr B2(M)(F1-F4). Five submillimeter components, 
Sgr B2(M)-SMA1 to SMA5, are detected in the central region 
of Sgr B2(M) in our observations. Outside the central region, 
seven components, Sgr B2(M)-SMA6 to SMAl 2 are identified. 
Given the positional accuracies of our observations and the ob- 
servations by de Pree, Goss & Gaume (1998), the projected po- 
sitions of Sgr B2(M)-SMA2, SMA6, SMAll and SMAl 2 coin- 
cide with those of the UC HII regions Fl, F10.37, F10.30, and 
F10.318, respectively. No centimeter source is detected towards 
Sgr B2(M)-SMA10, which is located south-west of the central 
region and displays an extended structure. 

Multi-component Gaussian fits were carried out towards 
bofli flie Sgr B2(N) and Sgr B2 (M) clumps using the IMFIT 
task. The residual fluxes after fitting are 4 and 11 Jy for Sgr 
B2(N) and Sgr B2(M), respectively. The large residual error 
in Sgr B2(M) clump is most likely caused by its complicated 
source structure. The peak positions, deconvolved angular sizes 
(FWHM), peak intensities, and total flux densities of the con- 
tinuum components are summarized in Table 1. The total flux 
densities of the Sgr B2(N) and Sgr B2 (M) cores are ~58 and 
61 Jy, while the peak fluxes measured by the bolometer array 
LABOCA at the 12-m telescope APEX are 150 and 138 Jy in a 
18.2" beam, respectively (Schuller, priv. conrai.), indicating that 
~60% of the flux is filtered out and that our SMA observations 
only pick up the densest parts of the Sgr B2 cores. 

Assuming that the 850 jum continuum is due to optically thin 
dust emission and using an average grain radius of 0.1 )Um, grain 
density of 3 g cm"-', and a gas-to-dust ratio of 100 (Hildebrand 
1983; Lis, Carlstrom & Keene 1991), the mass and column den- 
sity can be calculated using the formulae given in Lis, Carlstrom 
& Keen (1991). We adopt Qiv) = 4x10'^ at 850 /im (Hildebrand 
1983; Lis, Carlstrom & Keene 1991) and a dust temperature of 
150 K (Carlstrom & Vogel 1989; Lis et al. 1993) in the calcula- 
tion. On the basis of flux densities at 1.3 cm and 7 mm (Gaume 
et al. 1995; de Pree, Goss & Gaume 1998; Rolffs et al. 2011 
), most K and F subcomponents (except for K2 and F3) at 7 
mm have fluxes less than or comparable with those at 22.4 GHz, 
which produces descending spectra and is indicative of optically 
thin Hn regions at short wavelengths. The contributions of the 
free-free emission to the flux densities of the submillimeter com- 
ponents are smaller than 0.7% for K2 and F3 and smaller than 

0. 1% for other components, and can safely be ignored. The es- 
timated clump masses and column densities are given in Table 

1. The flux density of the fourteen detected submilUmeter com- 
ponents ranges from 1.2 to 47 Jy, corresponding to gas masses 
from 71 to 2731 M©. The column densities are a few times 10^' 
cm"^. Under the Rayleigh-Jeans approximation, 1 Jy beam"' 
in our SMA observations corresponds to a brightness temper- 
ature of 1 13 K. The peak brightness temperatures of Sgr B2(M)- 
SMAl and Sgr B2(N)-SMA1 are 270 and 200 K respectively, 
which place a lower limit on the dust temperatures at the peak 
position of the continuum. In this paper, the column densities 
and masses are estimated by use of source-averaged continuum 
fluxes. Based on the model fitting (Lis et al. 1991, 1993), flie 
adopted Q(v), dust temperature, and optically thin approxima- 
tion are reasonable guesses for the Sgr B2(N) and Sgr B2(M) 
cloud complexes. The total masses of Sgr B2(N) and (M), de- 
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Fig. 1. Continuum maps of Sgr B2 at 850 yum, with a synthesized beam of 0"4x0''24, PA=14.4° (lower-right corner in each panel). 
The left panels shows the image of Sgr B2(N), with contour levels (-1, 1, . . . 14)x4cr (lcr=0.031 Jy beam"'). The cross symbols 
indicate the positions of UC Hii regions detected in 1.3 cm continuum (Gaume et al. 1995). The right panel present the image of 
Sgr B2(M), with contour levels (-1, 1, 2, 3, 4, 4.5, 5.5, 5, 6, 7, 8, 10, . . . 28)x4o- (la-=0.021 Jy beam"'). The cross symbols indicate 
the positions of UC Hii regions detected in 7 mm continuum (de Pree, Goss & Gaume 1998). In each panel, the filled circle symbols 
present the peak positions of the submillimeter continuum sources. 



Table 1. Properties of the Continuum Sources 



Source 


a (J2000.0) 


6 (J2000.0) 


Deconvolved size 


Peak Intensity 


Flux Density 














(Jy beam"') 


(Jy) 


(lO-Mo) 


(lO-^cm-2) 



SgrB2(N) 
SgrB2(N) 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M). 
Sgr B2(M)- 
Sgr B2(M)- 
Sgr B2(M) 



SMAl 
SMA2 
■SMAl 
'-SMA2 
-SMA3 
SMA4 
SMA5 
-SMA6 
-SMA? 
SMA8 
SMA9 
SMAIO 
SMAll 
SMAl 2 



17 47 19.889 
17 47 19.885 
17 47 20.157 
17 47 20.133 
17 47 20.098 
17 47 20.150 
17 47 20.205 
17 47 20.171 
17 47 20.11 
17 47 20.216 
17 47 20.238 
17 47 19.989 
17 47 20.106 
17 47 20.132 



-28 22 18.22 
-28 22 13.29 
-28 23 04.53 
-28 23 04.06 
-28 23 03.9 
-28 23 03.26 
-28 23 04.78 
-28 23 05.91 
-28 23 06.18 
-28 23 06.48 
-28 23 06.87 
-28 23 05.83 
-28 23 03.01 
-28 23 02.24 



i:'72 X i:'28(-7.7°) 

i:'44 X iy02(-24.5°) 
i:'51 xO:'59(10.6°) 
0:'98x0:'58(14.7°) 
0:'93 X 0:'3(-9.9°) 
0:'55x0:'32(-1.5°) 
0:'66x0:'52(15.9°) 

0;'75 X 0:'43(-24.3°) 
0;'72 x0:'37(14.9°) 

0^52 X 0y45(-63.4°) 
0;'83 x0:'42(23.9°) 
3:'06 X i:'0(-25.8°) 

0:'97 X 0:'48(-27.7°) 
0:'82x0:'52(0.9°) 



1.79±0.039 
0.601±0.031 
2.39+0.138 

1.75+0.12 
0.866+0.09 
0.483+0.021 
0.52+0.051 
0.71+0.067 

0.6+0.051 
0.31+0.048 
0.369+0.032 
0.183+0.015 
0.695+0.058 
0.394+0.018 



47.48+1.034 
10.12+0.521 
20.9+1.21 
12.51+0.859 
4.147+0.431 
1.443+0.063 
2.488+0.244 
3.05+0.288 
2.366+0.201 
1.246+0.193 
1.922+0.164 
4.321+0.363 
4.423+0.37 
2.629+0.121 



27.31+0.59 

5.82+ 0.38 

12.02+0.7 

7.19+0.49 

2.38+0.25 

0.83+0.04 

1.43+0.14 

1.75+0.17 

1.36+0.12 

0.71+0.11 

1.11+0.09 

2.49+0.29 

2.54+0.21 

1.51+0.08 



4.54+0.1 
1.45+0.1 
4.94+0.29 
4.63+0.32 
3.13+0.33 
1.73+0.08 
1.53+0.15 
1.99+0.19 
1.87+0.16 
1.12+0.17 
1.16+0.1 
0.3+0.03 
2+0.17 
1.3+0.07 



Note: Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 



termined by summing up over the components, are comparable, 
3313 and 3532 M0, respectively, in spite of very different mor- 
phologies. We consider the adopted gas temperature of 150 K 
a reasonable guess for the average core temperatures, although 
some of the massive cores show higher peak brightness temper- 
atures. For those, the optically thin assumption is probably also 
not justified, and we would be underestimating their masses. 

4. Modeling 

A striking feature in the maps is the appearance of Sgr B2(N)- 
SMAl, which, although well resolved, does not appear to be 
fragmented. We used the three-dimensional radiative-transfer 
code RADMC-3EB developed by C. Dullemond, to model the 

^ http://www.ita.uni-heidelberg.de/"'dullemond/software/radmc-3d 



continuum emission of Sgr B2(N)-SMA1. Figure 2 shows three 
example models, whose radial profiles were obtained by Fourier- 
transforming the computed dust continuum maps, 'observing' 
with the uv coverage of the data, and averaging the image in cir- 
cular annuli. All models are heated by the stars in the UC Hii 
regions K2 and K3, assumed to have luminosities of ^ 10^ Lq, 
which is uncertain, but should give the right order of magni- 
tude (Rolffs et al. 2011). The dust mass opacity (0.6 cm^ g"') 
is interpolated from Ossenkopf & Henning (1994) without grain 
mantles or coagulation, which corresponds to a Q(v) of a few 
times 10"^ at 850 fim and is consistent with the value used in 
our calculation of masses and column densities. The model with 
a density distribution that follows the Plummer profile given by 

n = 1.7 X 10^ X ^1 -H ( iisooAu )^) ^2 cm"^ (half-power radius 
6500 AU) provides the best fit. We slso show in Fig. 2 a Gaussian 
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model with central density 2 x 10^ H2 cm ^ and half -power ra- 
dius 6500 AU, and a model whose density follows a radial power 

law, n= 10'^ X ( loooAu ) (outside of 1000 AU, the ra- 

dius of the Hii region K2). The latter model reproduces the peak 
flux measured by the bolometer array LABOCA at the 12-m tele- 
scope APEX, which is 150 Jy in a 18.2" beam, but does not fit 
the inner regions observed by the SMA very well. The Plummer 
and Gauss models fit Sgr B2(N)-SMA1 well, in terms of both 
shape and absolute flux. This implies that an extended compo- 
nent exists, which is filtered out in the SMA map but picked up 
by LABOCA. 




Radius [arcsec] 



Fig. 2. Radial profile of the Sgr B2(N)-SMA1 component, with 
errorbars denoting the rms in circular annuli. Overlaid are three 
models with different density distributions (red: Plummer; blue: 
Gaussian; green: power-law). The beam is depicted as a dashed 
Gaussian. 

It remains unclear whether the Plummer profile, which is 
also used to describe density profiles of star clusters, has any 
relevance to an evolving star cluster, or whether the cluster loses 
its memory of the gas density profile in the subsequent dynami- 
cal evolution. The model has a mass of around 3000 M© inside a 
radius of 1 1500 AU (0.056 pc), which represents an average den- 
sity of 4 X 10^ M0 pc"^. This is only the mass contained in the 
gas, and does not include the mass of the already formed com- 
pact objects providing the luminosity. Sgr B2(N) might represent 
a very young, embedded stage in the formation of a massive star 
cluster. 

5. Discussion 

Our SMA observations resolve the Sgr B2(M) and (N) cores 
into fourteen compact submillimeter continuum components. 
The two cores display very different morphologies. The source 
Sgr B2(N)-SMA1 is located north-east of the centimeter source 
K2, with an off'set of ~0.2". The continuum observations at 
1.3 cm having a resolution of 0.25" (Gaume et al. 1995) and 
at 7 mm a resolution of 0.1" (Rolffs et al. 2011) also detected 
a compact component centered on K2. We failed to detect any 
submillimeter continuum emission associated with sources Kl, 
K3, and K4. Sgr B2(N)-SMA2 is a high-mass dust core. 

In contrast to Sgr B2(N), a very fragmented cluster of high 
mass submillimeter sources is detected in Sgr B2(M). In addition 



to the two brightest and most massive components, Sgr B2(M)- 
SMAl and Sgr B2(M)-SMA2, situated in the centi-al region of 
Sgr B2(M), ten additional sources are detected, which indicates 
that there has been a high degree of fragmentation. The sensi- 
tivity of 0.021 Jy beam"' in our observations corresponds to a 
detectable gas mass of 1.2 Mq, but the observations are likely 
to be dynamic -range-limited, so it is difficult to determine the 
clump mass function in Sgr B2(M) down to smaller masses. 

The estimated column densities (10^^cm"^=33.4 gcm"^) in 
both the homogeneous starforming region Sgr B2(N) and the 
clustered Sgr B2(M) region are well in excess of the threshold 
of 1 g cm"^ for preventing cloud fragmentation and formation 
of massive stars (Krumholz & McKee 2008). The source sizes 
and masses derived from Gaussian fitting, assuming a spheri- 
cal source, lead to volume densities in excess of lO^cm"^ for all 
submillimeter sources detected in the SMA images. Assuming 
a gas temperature of 150 K, the thermal Jeans masses are less 
than 10 Mq, but the turbulent support is considerable. The large 
column densities, the gas masses, and the velocity field (Rolff's 
et al. 2010) suggest that the submillimeter components in the 
two regions are gravitationally unstable and in the process of 
forming massive stars. This process seems more advanced in 
Sgr B2(M), which is also reflected in the large number of em- 
bedded UC Hii regions found there. However, star formation in 
Sgr B2(N) does not appears to have progressed very far, a con- 
clusion also supported by the presence of only one UC Hii region 
embedded in one of the two clumps studied here. The observa- 
tions also showed that massive star formation taking place in the 
two clumps with outflow ages of ~ lO^* and ~ 10"* years for 
Sgr B2(N) and Sgr B2(M) clumps, respectively (Lis et al. 1993). 
If one can generalize these two examples, and if they provide 
snapshots in time of the evolution of basically equal cores, it 
seems that a massive star forms first in a relatively homogeneous 
core, another example with Plummer density profile and without 
fragmentation being the high-mass core G 10.47 (Rolffs et al., in 
prep.), followed by fragmentation or at least visible break-up of 
the core and subsequent star formation, perhaps aided by radia- 
tive or outflow feedback from the first star. This scenario may 
well apply to extremely high-mass cluster forming cores or to 
special environments only, since high-mass IRDCs have been 
shown to fragment early (Zhang et al. 2009). 
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